Open Access 



Research 



BMJ Open Postnatal epigenetic modification of 

glucocorticoid receptor gene in preterm 
infants: a prospective cohort study 

Masato Kantake, 1 Hiroshi Yoshitake, 2 Hitoshi Ishikawa, 3 Yoshihiko Araki, 2 
Toshiaki Shimizu 4 



To cite: Kantake M, 
Yoshitake H, Ishikawa H, 
etal. Postnatal epigenetic 
modification of glucocorticoid 
receptor gene in preterm 
infants: a prospective cohort 
study. BMJ Open 201 4;4: 
e005318. doi:10.1136/ 
bmjopen-201 4-00531 8 

► Prepublication history and 
additional material is 
available. To view please visit 
the journal (http://dx.doi.org/ 
10.1136/bmjopen-2014- 
005318). 



Received 22 March 2014 
Revised 5 June 2014 
Accepted 25 June 2014 



CrossMark 



For numbered affiliations see 
end of article. 



ABSTRACT 

Objective: To examine the environmental effects on 
cytosine methylation of preterm infant's DNA, because 
early life experiences are considered to influence the 
physiological and mental health of an individual 
through epigenetic modification of DNA. 
Design: A prospective cohort study, comparison of 
epigenetic differences in the glucocorticoid receptor 
(GR) gene between healthy term and preterm infants. 
Setting: Neonatal Intensive Care Unit in a Japanese 
University Hospital. 

Participants: A cohort of 40 (20 term and 20 
preterm) infants was recruited on the day of birth, and 
peripheral blood was obtained from each infant at birth 
and on postnatal day 4. 

Main outcome measures: The methylation rates in 
the 1-F promoter region of the GR gene using the 
Mquant method. 

Results: The methylation rate increased significantly 
between postnatal days 0 and 4 in preterm infants but 
remained stable in term infants. Thus, the methylation 
rate was significantly higher in preterm than in term 
infants at postnatal day 4. Several perinatal parameters 
were significantly correlated with this change in the 
methylation rate. Logistic regression analysis revealed 
that methylation rates at postnatal day 4 predicted the 
occurrence of later complications that required 
glucocorticoid administration during the neonatal 
period. No gene polymorphism was detected within the 
GR promoter region analysed. 
Conclusions: Although further large-scale studies are 
needed to detect the environmental factors that explain 
the difference in epigenetic modification among infants 
after birth, our data show that the postnatal 
environment influences epigenetic programming of GR 
expression through methylation of the GR gene 
promoter in premature infants, which may result in 
relative glucocorticoid insufficiency during the 
postnatal period. 
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INTRODUCTION 

Many human and animal studies have docu- 
mented the impact of early life experiences 
on the neurobiological regulation of stress 



Strengths and limitations of this study 



Postnatal environment resets the epigenetic 
modification of GR gene promoter during four 
days after birth in preterm infants. 
This epigenetic modification may result in relative 
glucocorticoid insufficiency during the postnatal 
period. 

Large-scale studies will be expected to clarify the 
environmental factors that explain the difference 
in epigenetic modification among infants after 
birth. 



responsiveness, mood and anxiety disor- 
ders. 1-4 Maternal care, familial functioning 
and childhood adversity all contribute to 
neurobiological regulation through epigen- 
etic modification of DNA, which is thought to 
be highly stable across a lifespan. 5-7 
Hippocampal glucocorticoid receptors (GRs) 
play a central role in modulating hypothal- 
amus pituitary adrenal (HPA) axis activity. 8-10 
Exon 1-F of the human GR gene is similar to 
rat exons 1-7, which is related to a maternal 
effect on cytosine methylation and expres- 
sion. 11 Several human studies have reported a 
relationship between early life experience 
and epigenetic modification of the GR 1-F 
promoter region. 11-16 These studies have 
shown that adverse early life experiences, 
such as maternal depression in the third tri- 
mester, disruption of adequate nurturing and 
history of childhood abuse, are associated 
with increased methylation of the GR pro- 
moter in leucocytes, 12-15 the hippocampus 11 
and the placenta. 16 

The stress of birth exceeds that of any 
other critical life event. Indeed, exposure to 
preterm labour and postnatal environmental 
stress, including the requisite acute care and 
prolonged physical separation under neo- 
natal intensive care, may be the most adverse 
environments during early life. Preterm birth 
is associated with significant long-term neu 
rodevelopmental impairments. 
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Increased methylation of the GR promoter in leuco- 
cytes may lead to down regulation of cell surface GR 
expression that, in turn, results in GR resistance (GCR). 
GCR refers to a decrease in the sensitivity of immune 
cells to glucocorticoid hormones that normally termin- 
ate the inflammatory response. 22-24 Evidence for GCR 
in response to prolonged stress has been found in 
humans, such as families of patients with cancer or 
those reporting high levels of loneliness. 25-27 On the 
other hand, relative adrenal insufficiency occurs when 
the HPA axis produces insufficient Cortisol relative to 
the degree of illness or stress. Several studies have 
reported a relative adrenal insufficiency in infants who 
develop chronic lung disease (CLD) 28 29 and cardiovas- 
cular instability 30-32 including late onset circulatory col- 
lapse (LCC). LCC is thought to be caused by late-onset 
adrenal insufficiency. An endocrinological study noted 
that even though relative glucocorticoid insufficiency in 
preterm infants is an underlying factor contributing to 
LCC, Cortisol concentrations did not differ between the 
LCC and control infant groups. 33 34 Therefore, another 
factor(s), including GCR, may be the cause of LCC. 

To test the hypothesis that an increase in cytosine 
methylation within the GR promoter during the peri- 
natal period in preterm infants causes physiological 
impairment in later life, we first compared epigenetic 
differences in the leucocyte GR promoter 1-F region 
(containing 33 CpG sites ranging from —3479 bp to 
—3259 bp upstream of the ATG start site) between 
healthy term and preterm infants on day 0 of birth and 
postnatal day 4. Furthermore, we examined the relation- 
ship between the postnatal change in the methylation 
rates and perinatal parameters and whether the methyla- 
tion status on postnatal day 4 caused complications in 
the neonatal period, such as CLD and LCC, which are 
thought to be caused by relative adrenal insufficiency. 

METHODS 

Prospective cohort study 

We recruited a cohort of 40 consecutive infants (20 term 
and 20 preterm) on the day of birth between April and 
October 2012 at our hospital. Eligible term infants were 
born within 37-41 weeks of gestation and received 
routine postnatal care as follows: day 0, skin to skin 
contact for 30 min from 30 min after birth; day 1, 
nursing four times (vaginal delivery (VD)) or one time 
(Caesarean section (C/S)); day 2, rooming with the 
baby (VD), nursing one time (C/S); day 3, rooming 
with the baby (VD), nursing six times (C/S); day 4, dis- 
charge (VD), rooming with the baby (C/S); day 5, 
rooming with the baby (C/S); day 6, discharge (C/S). 

An infant was defined as preterm if born before week 
37 of gestation and admitted to the neonatal intensive 
care unit (NICU) as premature. Patients born before 
37 weeks of gestation and who did not require admission 
to the NICU were excluded. All preterm infants in our 
study were exposed to maternal separation for at least 



4 days after birth. Small for gestational age (SGA) was 
defined as weight below the 10th centile for gestational 
age. The criteria for CLD were the need for additional 
oxygen after 28 days of age. 35 The criteria for LCC were 
(1) sudden systemic hypotension (mean arterial pressure 
<35 mm Hg sustained for at least 3h); (2) oliguria 
(<1 mL/kg/ h during an 8 h interval); (3) electrolyte 
abnormality (hyponatremia <130 mEq/L, hyperkalemia 
>5.5 mEq/L) and (4) excess body weight gain per day 
(15 g/kg per day or 1.5% per day). When more than 
two of these symptoms were detected seven or more 
days postnatally, and after all other common conditions 
causing systemic hypotension and oliguria, such as hypo- 
volemia, dehydration, symptomatic patent ductus arterio- 
sus and sepsis were ruled out, corticosteroid treatment 
was started. LCC was diagnosed once the infant 
responded to corticosteroid therapy. Daily changes in 
body weight, mean arterial pressure, serum sodium and 
potassium concentrations and urine output before and 
at symptom onset were also assessed to help the attend- 
ing physician make the diagnosis of LCC. 34 

Methylation analysis 

Peripheral blood was obtained from the umbilical cord 
(day 0) or peripheral vein (day 4) from each infant and 
stored at — 80° C until use. Blood was sampled at day 4 in 
term infants in conjunction with a metabolic screening 
test recommended for all infants born in Japan. 
According to the manufacturer's instruction, genomic 
DNA was extracted from 2 uL of frozen peripheral blood 
and bisulfite treated using an EZ DNA Methylation 
Direct Kit (Zymo Research Corp., Irvine, California, 
USA); 1 uL of the resulting 10 uL bisulfite-treated 
genomic DNA solution was then subjected to PCR to 
amplify the GR promoter 1-F region, as described previ- 
ously. 11 The Mquant method described by Leakey et aP 6 
was used to calculate the GR promoter methylation rate 
(see online supplementary methods). 

Statistical analysis 

We performed Wilcoxon's signed-rank test to evaluate 
the longitudinal difference between days 0 and 4 after 
birth in preterm and term group infants. We performed 
the Mann-Whitney U test to compare the methylation 
rates in the preterm and term groups at each CpG site. 
We performed the % 2 tests (for gestational age, birth 
weight and Apgar score (at 1 and 5 min)), and 
Mann-Whitney U tests (for intrauterine growth retard- 
ation, C/S delivery, respiratory distress, mechanical venti- 
lation, intracranial haemorrhage, bacterial infection and 
antenatal steroid administration) to compare the 
characteristics of each group of infants as shown in 
table 1. Stepwise multiple regression analysis was per- 
formed to investigate the relationships of the increased 
ratio of methylation rates as a dependent variable, with 
the perinatal parameters in preterm infants as inde- 
pendent variables. To perform this stepwise multiple 
regression analysis, several nominal variables were 
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Tablel Participant characteristics 
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p value 
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Gestational age (weeks) (mean±SD) 


30.8±3.2 


39.8±1.3 


<0.001 


Birth weight (g) (mean±SD) 


1431 ±605 


3076±271 


<0.001 


Intrauterine growth retardation (n) 


7 


1 


0.014 


Caesarean section delivery (n) 


15 


3 


<0.001 


Apgar score at 1 min (mean±SD) 


6.3±2.6 


9±0 


<0.001 


Apgar score at 5 min (mean±SD) 


8.4±1.2 


10±0 


<0.001 


Respiratory distress (n) 


10 


0 


<0.001 


Mechanical ventilation (n) 


10 


0 


<0.001 


Intracranial haemorrhage (n) 


0 


0 


1 


Bacterial infection (n) 


0 


0 


1 


Antenatal steroid administration (n) 


10 


0 


<0.001 



converted to categorical variables by grouping values 
into two categories. The adjusted coefficient of deter- 
mination (R 2 ) is the fraction of information of the 
dependent variable that is explained by the independent 
variables. Logistic regression analysis was performed to 
investigate the relationships between the methylation 
rate at 4 days as a dependent variable with later compli- 
cations in the neonatal period as independent variables. 
The goodness of fit of the logistic regression models was 
assessed using Pearson residuals. The pseudo-R 2 is the 
fraction of information of the dependent variable that is 
explained by the independent variables. 37 We used the 
likelihood ratio test to test the hypothesis regarding 
which dependent variable is significantly better among 
the four CpGs. A p value <0.05 (two-tailed) was consid- 
ered significant. All statistical analyses were carried out 
using the StatView R 5.0J for Windows (SAS Institute, 
Cary, North Carolina, USA) software package. 

RESULTS 

Participant characteristics 

Thirty-nine of the 40 participants completed the ana- 
lysis. Bisulflte-treated genomic DNA from one infant 
from the preterm group could not be amplified by PCR 
for unknown reasons. The characteristics of the infants 
are shown in table 1. Ten of 19 preterm infants received 
antenatal steroids. The preterm infants had higher rates 
of intrauterine growth retardation than that of infants 
delivered by G/ S. They had lower Apgar scores at 1 and 
5 min compared with those of term infants. Preterm 
infants displayed higher rates of respiratory distress and 
required mechanical ventilation more frequently during 
the postnatal period than did term infants. None of the 
infants experienced a bacterial infection or intracranial 
haemorrhage through postnatal day 4. No gene poly- 
morphism was detected within the GR promoter region 
analysed. 

Longitudinal changes in methylation status 

We examined 33 CpG sites in the GR 1-F promoter. In 
preterm infants, methylation rates significantly increased 



between birth and postnatal day 4 at 11 CpG sites 
(1, p=0.022; 2, p=0.044; 8, p=0.00029; 9, p=0.049; 10, 
p=0.044; 14, p=0.0013; 16, p=0.0022; 25, p=0.036; 26, 
p=0.049; 28, p=0.00097; 29, p=0.016) and significantly 
decreased at one CpG site (4, p=0.030). In contrast, 
methylation rates between birth and postnatal day 4 
were stable in term infants ( figure 1 ) . 

Cross-sectional differences in methylation status 

Methylation rates in preterm infants at birth were lower 
at three CpG sites (1; p=0.033, 5; p=0.0024, 8; p=0.0001) 
and higher at one site (4; p=0.040) compared with those 
in term infants. Methylation rates on postnatal day 4 
were significantly higher in preterm infants compared 
with term infants at seven CpG sites (15, p=0. 00061; 16, 
p=0.0083; 21, p=0.043; 25, p=0.0098; 26, p=0.0050; 27, 
p=0.00040; 28, p=0.038). Methylation rates in preterm 
infants were not lower than those of term infants at any 
CpG site (figure 2). 

Figure 3 shows the changes in methylation status in 
each group between birth and postnatal day 4 at CpG 
sites 16, 25, 26 and 28. Methylation rates at these sites 
significantly increased only in preterm infants, resulting 
in a significantly higher methylation rate in preterm 
than in term infants on postnatal day 4. 

Relationship between methylation rates at birth 
and perinatal parameters 

We analysed the relationship between the methylation 
rates at birth at CpG sites 1, 4, 5 and 8 and prenatal 
parameters such as gestational age, birth weight, SGA, 
antenatal steroid administration and mode of delivery. 
We found no significant correlation among them. 

Relationship between the increase in the methylation 
rate ratio and perinatal parameters 

To examine the relationship between the increased methy- 
lation rate ratio and perinatal parameters, we chose 12 
CpG sites (1, 2, 4, 8, 9, 10, 14, 16, 25, 26, 28, and 29) in 
which the methylation rates were significantly increased 
(with the exception of CpG4), and calculated the ratio as 
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the methylation rate at day 4 divided by that at day 0. Then 
we investigated the relationship between the change in 
ratio at each CpG site and the perinatal parameters of ges- 
tational age, birth weight, Apgar score, mode of delivery, 
SGA, respiratory distress syndrome (RDS), mechanical 
ventilation, antenatal steroid administration and admission 
to the NICU. We found associations between these para- 
meters and an increase in the methylation rate ratio at 
CpG4, CpG8, CpGlO, CpG25 and CpG29. The change in 
methylation ratio at each CpG site was associated with 
following perinatal parameters in the parenthesis: CpG4 
(gestational age and Apgar score at 1 min (regression 
coefhcient=0.029 and —0.085, respectively; overall 
p value=0.004; adjusted-R 2 =0.223) ) ; CpG8 (SGA, C/S, 
Apgar score at 5 min, and NICU admission (regression 
coefficient=0.43, —0.31, 0.16, and 0.79, respectively; overall 
p value <0.0001; adjusted-R 2 =0.573)); CpGlO (NICU admis- 
sion (regression coefhcient=0.29; overall p value=0.037; 
adjusted-R 2 =0.089)); CpG25 (SGA (regression coeffi- 
cient=0.86, overall p value=0.009, adjusted-R 2 =0.15)) and 
CpG29 (SGA (regression coefhcient=0.096, overall 
p value=0.039, adjusted-R 2 =0.086)), respectively. 



Relationship between methylation rate at postnatal day 4 
and complications during the neonatal period 

We chose four CpG sites (16, 25, 26 and 28) shown in 
figure 3 to investigate the relationship between methyla- 
tion status on postnatal day 4 and complications during 
the perinatal period. Then, we analysed the relationship 
between the methylation rate at these CpG sites on post- 
natal day 4 and subsequent complications. There were 
four preterm infants who had some complications and 
needed glucocorticoid administration during the neonatal 
period. (CLD; n=3, LCC; n=3) Except for these four cases, 
no infant received glucocorticoid therapy postnatally. 

The logistic regression models were statistically fit 
(Pearson residuals, x 2 =19.284, p=0.9799, pseudo-R 2 =0.429) . 
CpG16 had a significantly better explanatory factor among 
four CpGs (likelihood ratio test % 2 =3.889, p=0.0486) . 

DISCUSSION 

Our study is the first longitudinal analysis of GR pro- 
moter methylation in humans during the perinatal 
period, the stage that is thought to be the most 
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Figure 1 Longitudinal analysis of the methylation rate of the glucocorticoid receptor (GR) gene promoter in the leucocyte 
calculated using the Mquant method. The percentage of methylation at each CpG site (mean±SEM) in the GR promoter isolated 
in peripheral blood collected from term and preterm infants at birth (dark grey bar) and postnatal day 4 (light grey bar), 
respectively. Wilcoxon's signed-rank test was performed to compare the methylation rate on day 0 and day 4 at each CpG site. 
A p value <0.05 (two tailed) was considered statistically significant. Closed star symbol, significant increase in methylation rate; 
Open star symbol, significant decrease in methylation rate. No significant changes in methylation rate between postnatal days 0 
and 4 were observed in term infants. 
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Figure 2 Cross-sectional analysis of the methylation rates of the glucocorticoid receptor (GR) gene promoter in leucocytes 
determined using the Mquant method. Percentage of methylation at each CpG site (mean±SEM) in the GR promoter isolated 
from peripheral blood collected at birth and on postnatal day 4 in term (dark grey bar) and preterm (light grey bar) infants, 
respectively. Mann-Whitney U tests were performed to compare the methylation rates in preterm and term infants at each CpG 
site. A p value <0.05 (two tailed) was considered statistically significant. Closed star symbol, methylation rate significantly higher 
in preterm compared with term infants; Open star symbol, methylation rate significantly lower in preterm compared with term 
infants. 



susceptible to epigenetic regulation. The present study 
demonstrated that postnatal exposure to an adverse 
environment in preterm infants influences epigenetic 
programming of the glucocorticoid response via GR pro- 
moter methylation during the neonatal period. 

We found that perinatal parameters such as gestational 
age, intrauterine growth, Apgar score and mode of deliv- 
ery were predictors for changes in methylation rates at 
CpG sites in the GR promoter as well as NICU admission. 

It should be noted that gestational age was positively 
associated with an increase in the methylation ratio at 
CpG4. Exposure to preterm labour may be an adverse 
environment in early life. Therefore, the methylation 
rate ratio at CpG4 should be negatively associated with 
gestational age. CpG4 was the only site at which the 
methylation rates in preterm infants were significantly 
higher than those in term infants at birth. Moreover, 
CpG4 was the only site at which the methylation rate 
was lower on postnatal day 4 than on day 0. Taken 
together, the change in methylation at CpG4 seemed to 
be regulated inversely, compared with other CpG 
sites. The mechanisms underlying DNA methylation and 
demethylation probably develop due to unfavourable 



experiences during early life; however, recent studies 
that have investigated other transcripts of the GR gene 
have proposed a more complex child maltreatment- 
induced DNA methylation model than previously 
described for exon 1-F. 38 In the present study, CpG4 may 
be a transcriptionally silent methylation site that 
mediated aberrant epigenetic regulation, as proposed by 
Lutz et al. 38 Low Apgar score, which is thought to be an 
unfavourable experience during the perinatal period, is 
positively associated with the CpG4 methylation ratio 
increase and negatively associated with that of CpG8. 
Thus, low Apgar score may have a negative effect on 
methylation in most CpG sites in the GR 1-F promoter. 
One possible explanation is that low Apgar score indi- 
cates an adverse condition of cardiopulmonary function, 
which, in turn, inactivates the enzymatic cascade for 
DNA methylation. Stable oxygen and energy supple- 
ments may be necessary to open the window for epigen- 
etic programming by DNA methylation. C/S delivery was 
negatively associated with an increase in the methylation 
ratio at CpG8. This obstetrical intervention may partially 
reduce stress-induced epigenetic programming of the 
GR promoter. SGA was positively associated with an 
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Figure 3 Change in methylation status between birth and postnatal day 4 at CpG sites 16, 25, 26 and 28. Methylation rates at 
these sites were significantly greater in preterm than in term infants resulting in a significantly higher methylation rate ratio in the 
preterm infants on postnatal day 4. The day 0 and day 4 methylation rate ratios were compared using Wilcoxon's signed-rank 
test. The Mann -Whitney U test was performed to compare the methylation rates in preterm infants with those in term infants, ns, 
not significant. 



increase in methylation rate ratio at CpG8, 25 and 29. 
The intrauterine environment, like the postnatal envir- 
onment, can affect fetal epigenetic modification of 
DNA. SGA infants have been thought to demonstrate 
epigenomic dys regulation, which mediates the long-term 
consequences of intrauterine growth restriction at 
birth. 39 40 The present study shows that, in the GR gene, 
SGA might influence the infant's susceptibility to the 
postnatal environment rather than epigenomic dysregu- 
lation at birth. As we expected, NICU admission fol- 
lowed by maternal separation was positively associated 
with an increase in the methylation rate ratio at CpG8 
and 10. Taken together, we could approximately estimate 
the relationship between perinatal parameters and 
changes in methylation rates at CpG sites in the GR pro- 
moter. Further precise analysis should give insights into 
the host-environment interaction. 

Four preterm infants needed glucocorticoid adminis- 
tration during the postnatal period due to CLD or LCC. 
The logistic regression models were statistically fit to 
predict the development of these complications (CLD 
and LCC), which were thought to be partly related to 



relative adrenal insufficiency. CpG16 had a significantly 
better explanatory factor among the four CpGs in 
which methylation rates increased significantly only in 
preterm infants resulting in a significantly higher methy- 
lation rate in preterm than in term infants on postnatal 
day 4. CpG16 is a binding site for nerve growth factor 
inducible protein A (NGFI-A), which has emerged as a 
central regulator of early inflammatory and immune pro- 
cesses and potentiates GR 1-F promoter activity. 11 41-43 

Taken together, our results suggest that the high methyla- 
tion rate observed in leucocyte GR promoters may result in 
GR unresponsiveness to early inflammatory stimuli by 
NGFI-A, followed by GCR and/ or a failure to downregulate 
inflammatory responses by glucocorticoids. 44 

One limitation of our study was the small sample size. 
Further studies with larger sample sizes that include 
patients with chronic inflammation are needed. 



CONCLUSIONS 

We found marked changes in DNA methylation induced 
by the postnatal environment during the 4 days 
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immediately after birth. Moreover, the methylation status 
at postnatal day 4 may predict the development of com- 
plications during the neonatal period, which is thought 
to be induced by relative glucocorticoid insufficiency in 
preterm infants. 
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